The EH domain-containing protein, EdeA, is involved in endocytosis, cell wall integrity, and pathogenicity in Aspergillus fumigatus.

Endocytosis has been extensively studied in yeasts, where it plays crucial roles in growth, signaling regulation, and cell-surface receptor internalization. However, the biological functions of endocytosis in pathogenic filamentous fungi remain largely unexplored. In this study, we aimed to functionally characterize the roles of EdeA, an ortholog of the Saccharomyces cerevisiae endocytic protein Ede1, in Aspergillus fumigatus. EdeA was observed to be distributed as patches on the plasma membrane and concentrated in the subapical collar of hyphae, a localization characteristic of endocytic proteins. Loss of edeA caused defective hyphal polarity, reduced conidial production, and fewer sites of endocytosis initiations than that of the parental wild type. Notably, the edeA null mutant exhibited increased sensitivity to cell wall-disrupting agents, indicating a role for EdeA in maintaining cell wall integrity in A. fumigatus. This observation was further supported by the evidence showing that the thickness of the cell wall in the ΔedeA mutant increased, accompanied by abnormal activation of MpkA, a key component in the cell wall integrity pathway. Additionally, the ΔedeA mutant displayed increased pathogenicity in the Galleria mellonella wax moth infection model, possibly due to alterations in cell wall morphology. Site-directed mutagenesis identified the conserved residue E348 within the third EH (Eps15 homology) domain of EdeA as crucial for its subcellular localization and functions. In conclusion, our results highlight the involvement of EdeA in endocytosis, hyphal polarity, cell wall integrity, and pathogenicity in A. fumigatus. IMPORTANCE: Aspergillus fumigatus is a significant human pathogenic fungus known to cause invasive aspergillosis, a disease with a high mortality rate. Understanding the basic principles of A. fumigatus pathogenicity is crucial for developing effective strategies against this pathogen. Previous research has underscored the importance of endocytosis in the infection capacity of pathogenic yeasts; however, its biological function in pathogenic mold remains largely unexplored. Our characterization of EdeA in A. fumigatus sheds light on the role of endocytosis in the development, stress response, and pathogenicity of pathogenic molds. These findings suggest that the components of the endocytosis process may serve as potential targets for antifungal therapy.

Electroacupuncture Improves Insulin Resistance in Type 2 Diabetes Mice by Regulating Intestinal Flora and Bile Acid.

Introduction: Adjusting internal organs and dredging channel electroacupuncture has a definite effect on type 2 diabetes, but the specific mechanism still needs to be further clarified. This study aims to investigate the effects of electroacupuncture on the gut microbiota and bile acids in db/db mice after the intervention of "adjusting internal organs and dredging channel" and further explore its mechanism of action in treating T2DM. Methods: We used db/db mice as the animal model and db/m mice from the same litter as the blank control group, a total of 4 weeks of intervention were conducted. We evaluated the effectiveness of the "adjusting internal organs and dredging channel" treatment by detecting indicators related to glucose and lipid- metabolism. Detect changes in the gut microbiota of mice in each group using 16SrDNA sequencing technology. The content of bile acids in mouse feces was determined using liquid chromatography mass spectrometry, and the correlation analysis between different bile acids and differential bacterial communities was performed. The expression levels of TGR5 and GLP-1 proteins were measured using the Western blot method. Results: Adjusting internal organs and dredging channel electroacupuncture can improve blood glucose levels in db/db mice, increase the abundance of Firmicutes and Actinobacteria, and increase the content of fecal bile acid pool heavy CA and UDCA. At the same time, it also increased the content of TGR5/GLP1 in the small intestine. Conclusion: Adjusting internal organs and dredging channel electroacupuncture can improve the disorder of glucose and lipid metabolism in db/db mice, regulate the abundance and colony composition of intestinal microbiota in mice, and regulate bile acid metabolism in mice. The interaction between bile acid and intestinal microbiota can also be observed; Mutual influence may play a role in regulating blood sugar together.

Transcription factor CreA is involved in the inverse regulation of biofilm formation and asexual development through distinct pathways in Aspergillus fumigatus.

The exopolysaccharide galactosaminogalactan (GAG) contributes to biofilm formation and virulence in the pathogenic fungus Aspergillus fumigatus. Increasing evidence indicates that GAG production is inversely linked with asexual development. However, the mechanisms underlying this regulatory relationship are unclear. In this study, we found that the dysfunction of CreA, a conserved transcription factor involved in carbon catabolite repression in many fungal species, causes abnormal asexual development (conidiation) under liquid-submerged culture conditions specifically in the presence of glucose. The loss of creA decreased GAG production independent of carbon sources. Furthermore, CreA contributed to asexual development and GAG production via distinct pathways. CreA promoted A. fumigatus GAG production by positively regulating GAG biosynthetic genes (uge3 and agd3). CreA suppressed asexual development in glucose liquid-submerged culture conditions via central conidiation genes (brlA, abaA, and wetA) and their upstream activators (flbC and flbD). Restoration of brlA expression to the wild-type level by flbC or flbD deletion abolished the abnormal submerged conidiation in the creA null mutant but did not restore GAG production. The C-terminal region of CreA was crucial for the suppression of asexual development, and the repressive domain contributed to GAG production. Overall, CreA is involved in GAG production and asexual development in an inverse manner.

Transcription factors SltA and CrzA reversely regulate calcium homeostasis under calcium-limited conditions.

IMPORTANCE: Calcium ions are ubiquitous intracellular signaling molecules for many signaling pathways regulating the fungal response to stress and antifungal drugs. The concentration of intracellular calcium is tightly regulated in its storage, release, and distribution. CrzA is the best-studied transcription factor that regulates this process under sufficient calcium or other external signals. However, CrzA was excluded from nuclei and then lost transcriptional activation under calcium-limited conditions. The regulators in the Ca2+ signaling pathway under calcium-limited conditions remain unclear. Here, we identified SltA as a key regulator in the Ca2+ signaling pathway under calcium-limited conditions, and the underlying mechanisms were further explored in Aspergillus fumigatus. These findings reveal a transcriptional control pathway that precisely regulates calcium homeostasis under calcium-limited conditions.

The C2H2-Type Transcription Factor ZfpA, Coordinately with CrzA, Affects Azole Susceptibility by Regulating the Multidrug Transporter Gene atrF in Aspergillus fumigatus.

The incidence of invasive aspergillosis caused by Aspergillus fumigatus has risen steadily over the past few decades due to the limited effective treatment options and the emergence of antifungal-resistant isolates. In clinic-isolated A. fumigatus, the azole resistance mechanism is primarily caused by mutations of the drug target and/or overexpression of drug efflux pumps. However, knowledge about how drug efflux pumps are transcriptionally regulated is limited. In this study, we found that loss of a C2H2 transcription factor ZfpA (zinc finger protein) results in the marked upregulation of a series of drug efflux pump-encoding genes, especially atrF, which contributes to azole drug resistance in A. fumigatus. CrzA is a previously identified positive transcription factor for genes of drug efflux pumps, and ZfpA transcriptionally inhibits expressions of drug efflux pumps in a CrzA-dependent way. Under the treatment of azoles, both ZfpA and CrzA transfer to nuclei and coregulate the expression of multidrug transporters and then keep normal drug susceptibility in fungal cells. Findings in this study demonstrated that ZfpA is not only involved in fungal growth and virulence potential but also negatively regulates antifungal drug susceptibility. IMPORTANCE Conserved across all kingdoms of life, ABC transporters comprise one of the largest protein families. They are associated with multidrug resistance, affecting aspects such as resistance to antimicrobials or anticancer drugs. Despite the importance of ABC transporters in multidrug resistance, the understanding of their regulatory network is still limited in A. fumigatus. Here, we found that the loss of the transcription factor ZfpA induces the expression of the ABC transporter gene atrF, altering azole susceptibility in A. fumigatus. ZfpA, coordinately with CrzA, affects the azole susceptibility by regulating the expression of the ABC transporter gene atrF. These findings reveal the regulatory mechanism of the ABC transporter gene atrF in A. fumigatus.

Case Report: Partial response to single-agent pembrolizumab in a chemotherapy-resistant metastatic pancreatic cancer patient with a high tumor mutation burden.

Single-agent immune checkpoint blockade has shown no clinical benefits in pancreatic cancer. Recently, the programmed cell death protein 1 (PD-1) antibody pembrolizumab has been recommended as a treatment option for high tumor mutational burden (TMB) solid tumors based on the data from a basket trial. However, no pancreatic cancer patients were enrolled in that trial. Whether pancreatic cancer patients with high TMB respond to PD-1 blockade as well remains unclear. Here, we report a case with a partial response to single-agent immunotherapy with pembrolizumab in pancreatic cancer with high TMB after the failure of several lines of chemotherapy. This result indicates that single-agent immunotherapy may be effective in pancreatic cancer patients with high TMB. In addition, in order to understand the basic immune state of our patients, we also analyzed the changes in immune cells in peripheral blood with cytometry by time-of-flight mass spectrometry (CyTOF) before and after pembrolizumab treatment.

Revealing diverse alternative splicing variants of the highly homologous SMN1 and SMN2 genes by targeted long-read sequencing.

The survival of motor neuron (SMN) genes, SMN1 and SMN2, are two highly homologous genes related to spinal muscular atrophy (SMA). Different patterns of alternative splicing have been observed in the SMN genes. In this study, the long-read sequencing technique for distinguishing SMN1 and SMN2 without any assembly were developed and applied to reveal multiple alternative splicing patterns and to comprehensively identify transcript variants of the SMN genes. In total, 36 types of transcript variants were identified, with an equal number of variants generated from both SMN1 and SMN2. Of these, 18 were novel SMN transcripts that have never been reported. The structures of SMN transcripts were revealed to be much more complicated and diverse than previously discovered. These novel transcripts were derived from diverse splicing events, including skipping of one or more exons, intron retention, and exon shortening or addition. SMN1 mainly produces FL-SMN1, SMN1Δ7, SMN1Δ5 and SMN1Δ3. The distribution of SMN2 transcripts was significantly different from those of SMN1, with the majority transcripts to be SMN2Δ7, followed by FL-SMN2, SMN2Δ3,5 and SMN2Δ5,7. Targeted long-read sequencing approach could accurately distinguish sequences of SMN1 from those of SMN2. Our study comprehensively addressed naturally occurring SMN1 and SMN2 transcript variants and splicing patterns in peripheral blood mononuclear cells (PBMCs). The novel transcripts identified in our study expanded knowledge of the diversity of transcript variants generated from the SMN genes and showed a much more comprehensive profile of the SMN splicing spectrum. Results in our study will provide valuable information for the study of low expression level of SMN proteins and SMA pathogenesis based on transcript levels.

The relationship between ecological executive function and stigma among patients with epilepsy: The mediating effect of social support.

OBJECTIVE: Stigma is a key determinant of mental health among patients with epilepsy (PWE). In prior work, many PWE have been found to exhibit impaired executive function (EF) and to lack sufficient social support. This study was developed to explore factors that may influence stigma in PWE with a focus on the associations among EF, social support, and stigma in this patient population. METHODS: A questionnaire was administered to 121 patients with primary epilepsy as a means of collecting clinical and demographic details. The severity of EF impairment, social support levels, and stigma were analyzed using the Behavior Rating Inventory of Executive Function-Adult (BRIEF-A), the Social Support Rating Scale, and the Kilifi Stigma Scale for Epilepsy-Chinese (KSSE-C). Descriptive analyses were used to evaluate demographic details, and parameters associated with stigma were identified through multiple linear stepwise analyses. Spearman's correlation analyses and moderated mediation analyses were utilized to examine relationships between impaired EF, social support, and stigma. Significant mediation effects were identified using the SPSS PROCESS macro via a bootstrap approach. RESULTS: Up to 95.9% of the 121 PWE in this study were affected by stigma, with differneces in stigma being associated with long-term residence (urban or rural), family monthly income, disease control status (controlled or uncontrolled), and medication types. Multiple stepwise linear regression analyses revealed that the number of years of education (t = - 4.58, P < 0.001), family monthly income per capita (t = -3.43, P = 0.001), and age of first episode (t = -2.71, P = 0.008) for PWE were negatively correlated with stigma scores, while long-term residence (t = 2.79, P = 0.006), course of disease (t = 3.65, P < 0.001), disease control (t = 2.79, P = 0.006), and types of medication (t = 2.73, P = 0.007) were positively correlated with stigma scores. Impaired EF was found to be significantly associated with stigma (ß = 0.64, P < 0.001), and social support was able to mediate the association between stigma and impaired EF (ß = 0.64, P < 0.001). Non-parametric bootstrap analysis results revealed the indirect impact of EF through social support (95% bootstrap CI = 0.11, 0.22), with the indirect influence of social support accounting for 44.3% of the overall effect of stigma. CONCLUSIONS: These results offer new insight into the interactive mechanisms that underlie stigma and EF, in addition to clarifying the ability of social support to mediate this relationship. These data thus highlight valuable theoretical and methodological approaches to preventing stigma in PWE, suggesting that stigma in PWE can be effectively reduced by bolstering EF, and social support plays a mediating effect in this context.

The sterol C-14 reductase Erg24 is responsible for ergosterol biosynthesis and ion homeostasis in Aspergillus fumigatus.

Ergosterol, a major lipid present in the fungal cell membrane, is considered as an effective antifungal drug target. A rational strategy for increasing drug reservoir relies on functionally validation of essential enzymes involved in fungal key biological pathway. Current knowledge regarding the essential genes in the ergosterol biosynthesis pathway is still limited in the opportunistic human pathogen Aspergillus fumigatus. In this study, we characterized two endoplasmic reticulum-localized sterol C-14 reductases encoded by both erg24A and erg24B homologs that are essential for the viability of A. fumigatus despite the fact that neither paralog is essential individually. Loss of one homolog of Erg24 impairs hyphal growth, conidiation, and virulence but has no effect on ergosterol biosynthesis. To investigate the functional significance of erg24, a conditional double mutant (Δerg24B niiA::erg24A) was constructed in the Δerg24B background. Strikingly, the conditional erg24 double mutant exhibited severe growth defects and accumulation of sterol intermediate. Moreover, the addition of metal ions and the overexpression of the corresponding ion transporters could rescue the growth defects of the erg24 double mutant in A. fumigatus, implying that the defective phenotype of the erg24 double mutant is tightly associated with dysregulation of ion homeostasis. Taken together, our results demonstrate the critical role of Erg24 in ergosterol biosynthesis and ion homeostasis in A. fumigatus, which may have important implications for antifungal discovery. KEY POINTS: • We characterized two endoplasmic reticulum-localized sterol C-14 reductases Erg24A and Erg24B in A. fumigatus. • Erg24A and Erg24B in combination, but not individually, are required for the viability of A. fumigatus. • Inactivation of Erg24 leads to the disruption of ion homeostasis and affects ergosterol biosynthesis.

Altered Inflammatory Pathway but Unaffected Liver Fibrosis in Mouse Models of Nonalcoholic Steatohepatitis Involving Interleukin-1 Receptor-Associated Kinase 1 Knockout.

BACKGROUND Interleukin-1 receptor-associated kinases (IRAKs) are crucial mediators in the signaling pathways of Toll-like receptors (TLRs)/IL1Rs. Targeting the IRAK4/IRAK1/TRAF6 axis and its associated pathway has therapeutic benefits in liver fibrosis. However, the function of IRAK1 itself in the development of liver fibrosis remains unknown. MATERIAL AND METHODS Irak1 global knockout (KO) mice were generated to study the functional role of Irak1 in liver fibrosis. Male Irak1 knockout and control mice were challenged with chronic carbon tetrachloride (CCl4) or fed a methionine- and choline-deficient diet (MCDD) to generate models of nonalcoholic steatohepatitis (NASH). Liver inflammation and collagen deposition were assessed by histological examination, quantitative real-time PCR (qRT-PCR), and western blotting of hepatic tissues. RESULTS The mRNA expression of the downstream inflammatory gene Il1b was significantly lower in Irak1-KO than in control mice. Irak1 ablation had little effect on inflammatory cell infiltration into livers of mice with NASH. Collagen deposition and the expression of genes related to fibrogenesis were similar in the livers of Irak1-KO and control mice exposed to CCl4 and MCDD. The loss of Irak1 did not affect lipid or glucose metabolism in these experimental models of steatohepatitis. CONCLUSIONS Irak1 knockout reduced the expression of inflammatory genes but had no effect on hepatic fibrogenesis. The Irak1-related pathway may regulate liver fibrosis via other pathways or be compensated for by other factors.

Biochemical and genetic approaches to the prenatal diagnosis of propionic acidemia in 78 pregnancies.

BACKGROUND: Propionic acidemia (PA) is a serious metabolic disorder, and different approaches have been applied to its prenatal diagnosis. To evaluate the reliability and validity of a biochemical strategy in the prenatal diagnosis of PA, we conducted a retrospective study of our 11-year experiences at a single center. METHODS: We accumulated data from 78 pregnancies from 58 families referred to our center and provided prenatal diagnosis by directed genetic analysis and/or metabolite measurement using tandem mass spectrometry (MS/MS) and gas chromatography/mass spectrometry (GC/MS) of amniotic fluid (AF) samples. RESULTS: Sixty-five unaffected fetuses (83.33%) and 13 affected fetuses (16.67%) were confirmed in our study. The characteristic metabolites including propionylcarnitine (C3) level, C3/acetylcarnitine (C2) ratio and 2-methylcitric acid (2MCA) level in unaffected and affected groups showed significant differences (P < 0.0001), while the level of 3-hydroxypropionic acid (3HPA) showed no significant difference between the two groups (P > 0.05).Of the 78 pregnancies, 24 fetuses were found to have either one causative pathogenic variant or were without genetic information in the proband. Three of these fetuses had elevated AF levels of C3, C3/C2 ratio, and 2MCA and, thus, were determined to be affected, while the remaining fetuses were determined to be unaffected based on a normal AF metabolite profile. Our genetic and biochemical results were highly consistent with postnatal follow-up results on all unaffected fetuses. CONCLUSIONS: We conclude that a biochemical approach can serve as a fast and convenient prenatal diagnostic method for pregnancies at an increased risk for PA, which could be used in conjunction with genetic testing for precise prenatal diagnosis of this disorder. In our analysis, the characteristic metabolites C3 level, C3/C2 ratio, and 2MCA level in AF supernatant were dependable biochemical markers for diagnosis, of which the C3/C2 ratio appears to be the most reliable biochemical marker for the prenatal diagnosis of PA.

Transcriptome analysis of Cryptococcus humicola under aluminum stress revealed the potential role of the cell wall in aluminum tolerance.

Aluminum (Al) toxicity is one of the most important limiting factors for crop yield in acidic soils. Bound Al gets converted into a toxic ionic state (Al3+) in acidic soil. Recent studies have shown that Al can act on the cell walls, cell membranes, organelles, and nuclei of microorganisms and affect substance and energy metabolism. To explore the gene expression at the transcriptional level under Al stress, we sequenced the transcriptome of Cryptococcus humicola, which is a highly Al-resistant yeast strain isolated from acidic soil and tolerates up to 200 mM Al3+. The screening conditions for genes from the control and experimental group were a false discovery rate (FDR) <0.05 and log 2|FC| > 1. A total of 4760 genes were differentially expressed, among which 3066 were upregulated and 1694 were downregulated. These genes control glycometabolism, protein synthesis, lipid metabolism and signalling pathways. Eleven selected differentially expressed genes were further validated using qRT-PCR. The results suggested that Al stress leads to complex responses in C. humicola. The effects of Al on the ß-d-glucan and mannose contents and Al accumulation in the cell wall were determined. With an increase in the Al treatment time and concentration, the contents of ß-d-glucan and mannose showed a trend of first increasing and then decreasing. Under Al treatment, the Al content of the cell wall also showed a trend of first increasing and then decreasing. These results suggested that Al accumulates in the cell wall and the cell wall plays a vital role in the Al resistance of C. humicola. The differentially expressed genes provide a foundation for the further study of Al tolerance in C. humicola.

Identification of novel SMN1 subtle mutations using an allelic-specific RT-PCR.

Spinal muscular atrophy (SMA) is caused by homozygous deletions of the SMN1 gene in approximately 95% of patients. The remaining 5% of patients with SMA retain at least one copy of the SMN1 gene carrying insertions, deletions, or point mutations. Although molecular genetic testing for most SMA patients is quite easy, diagnosing "nondeletion" SMA patients is still compromised by the presence of a highly homologous SMN2 gene. In this study, we analyzed the SMN1/SMN2 copy number by quantitative PCR and multiplex ligation-dependent probe amplification (MLPA). Further, common primers for both SMN1 and SMN2 sequences were used to screen DNA intragenic mutations. To confirm whether the identified mutations occurred in SMN1 or SMN2, we improved the traditional RT-PCR method by only amplifying SMN1 transcripts using an allelic-specific PCR (AS-RT-PCR) strategy. We identified six SMN1 point mutations and small indels in 8 families, which included c.683T>A, c.22dupA, c.815A>G, c.19delG, c.551_552insA and c.401_402delAG. To the best of our knowledge, the latter three have never been previously reported. The most common mutation in Chinese patients is c.22dupA, which was identified in three families. In this work, we demonstrated AS-RT-PCR to be reliable for identifying SMN1 subtle mutations, especially the prevalent mutation c.22dupA in Chinese SMA patients. By reviewing published papers and summarizing reported SMN1 mutations, a distinct ethnic specificity was found in SMA patients from China. Our research extends the SMN1 mutation spectrum.

Prenatal diagnosis of methylmalonic aciduria from amniotic fluid using genetic and biochemical approaches.

OBJECTIVES: This study reported the clinical prenatal diagnosis experience of families affected by methylmalonic acidemia (MMA) evaluated at a single prenatal diagnosis center over 8 years, and the reliability of a biochemical approach for prenatal diagnosis was analyzed. METHODS: Prenatal diagnosis data for 187 MMA families referred to our center from 2009 to 2016 were reviewed retrospectively. The results of the genetic analysis and biochemical approach were compared. RESULTS: A total of 41 MMA-affected pregnancies (21%) were identified. The biochemical analysis could identify the true status of 99.5% of fetuses. The diagnostic sensitivities of the propionylcarnitine (C3) level, the C3 to acetylcarnitine (C2) ratio (C3/C2), the methylmalonic acid, and methylcitrate levels in the amniotic fluid were 95.1%, 100%, 100%, and 82.9%, respectively, and the specificities were 98.7%, 99.3%, 97.4%, and 96.7%, respectively. CONCLUSIONS: The biochemical analysis could be optionally used in the prenatal diagnosis of MMA, especially in cases where the genetic results are inconclusive. Among the four tested biochemical markers, C3/C2 appeared to be the most reliable.

[Diagnosis of two cases from one family with Joubert syndrome caused by novel mutations of TCTN1 gene by whole exome sequencing].

OBJECTIVE: To explore the pathogenesis of two fetuses from one family affected with Joubert syndrome (JS). METHODS: Whole exome sequencing was employed to screen potential mutations in both fetuses. Suspected mutations were verified by Sanger sequencing. Impact of intronic mutations on DNA transcription was validated by cDNA analysis. RESULTS: Two novel TCTN1 mutations, c.342-8A>G and c.1494+1G>A, were identified in exons 2 and 12, respectively.cDNA analysis confirmed the pathogenic nature of both mutations with interference of normal splicing resulting in production of truncated proteins. CONCLUSION: The genetic etiology of the family affected with JS has been identified.Above findings have enriched the mutation spectrum of TCTN1gene and facilitated understanding of the genotype-phenotype correlation of JS.

Differential expression of Meis2, Mab21l2 and Tbx3 during limb development associated with diversification of limb morphology in mammals.

Bats are the only mammals capable of self-powered flight using wings. Differing from mouse or human limbs, four elongated digits within a broad wing membrane support the bat wing, and the foot of the bat has evolved a long calcar that spread the interfemoral membrane. Our recent mRNA sequencing (mRNA-Seq) study found unique expression patterns for genes at the 5' end of the Hoxd gene cluster and for Tbx3 that are associated with digit elongation and wing membrane growth in bats. In this study, we focused on two additional genes, Meis2 and Mab21l2, identified from the mRNA-Seq data. Using whole-mount in situ hybridization (WISH) we validated the mRNA-Seq results for differences in the expression patterns of Meis2 and Mab21l2 between bat and mouse limbs, and further characterize the timing and location of the expression of these two genes. These analyses suggest that Meis2 may function in wing membrane growth and Mab21l2 may have a role in AP and DV axial patterning. In addition, we found that Tbx3 is uniquely expressed in the unique calcar structure found in the bat hindlimb, suggesting a role for this gene in calcar growth and elongation. Moreover, analysis of the coding sequences for Meis2, Mab21l2 and Tbx3 showed that Meis2 and Mab21l2 have high sequence identity, consistent with the functions of genes being conserved, but that Tbx3 showed accelerated evolution in bats. However, evidence for positive selection in Tbx3 was not found, which would suggest that the function of this gene has not been changed. Together, our findings support the hypothesis that the modulation of the spatiotemporal expression patterns of multiple functional conserved genes control limb morphology and drive morphological change in the diversification of mammalian limbs.

Unique expression patterns of multiple key genes associated with the evolution of mammalian flight.

Bats are the only mammals capable of true flight. Critical adaptations for flight include a pair of dramatically elongated hands with broad wing membranes. To study the molecular mechanisms of bat wing evolution, we perform genomewide mRNA sequencing and in situ hybridization for embryonic bat limbs. We identify seven key genes that display unique expression patterns in embryonic bat wings and feet, compared with mouse fore- and hindlimbs. The expression of all 5'HoxD genes (Hoxd9-13) and Tbx3, six known crucial transcription factors for limb and digit development, is extremely high and prolonged in the elongating wing area. The expression of Fam5c, a tumour suppressor, in bat limbs is bat-specific and significantly high in all short digit regions (the thumb and foot digits). These results suggest multiple genetic changes occurred independently during the evolution of bat wings to elongate the hand digits, promote membrane growth and keep other digits short. Our findings also indicate that the evolution of limb morphology depends on the complex integration of multiple gene regulatory networks and biological processes that control digit formation and identity, chondrogenesis, and interdigital regression or retention.